Continuous Quantum Phase Transition in a Kondo Lattice Model 
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We study the magnetic quantum phase transition in an anisotropic Kondo lattice model. The 
dynamical competition between the RKKY and Kondo interactions is treated using an extended 
dynamic mean field theory (EDMFT) appropriate for both the antiferromagnetic and paramag- 
netic phases. A quantum Monte Carlo approach is used, which is able to reach very low temper- 
atures, of the order of 1% of the bare Kondo scale. We find that the finite-temperature magnetic 
transition, which occurs for sufficiently large RKKY interactions, is first order. The extrapolated 
zero-temperature magnetic transition, on the other hand, is continuous and locally critical. 

PACS numbers: 71.10.Hf, 71.27.-|-a, 75.20.Hr, 71.28.-(-d 
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The observation of magnetic quantum critical points 
and associated non-Fermi liquid behavior has raised a 
renewed interest in Kondo lattice systems Q. An anti- 
ferromagnetic phase transition in dimensions higher than 
one was anticipated very early on 0, 0, • For a spin 
S = Kondo lattice system that is not too far away from 
half-filling, an antiferromagnetic metal phase is expected 
if the RKKY interactions dominate over the Kondo in- 
teractions, whereas a paramagnetic heavy fermion phase 
should occur in the opposite limit. Except for some 
special cases however, there has been a lack of dy- 
namical treatments of the competition between antifer- 
romagnetism and Kondo effects in lattice problems. A 
key issue that has been left open is whether the zero- 
tcmpcrature magnetic phase transition is first order or 
continuous; certainly quantum critical fiuctuations can 
be of any relevance only if the transition is cither contin- 
uous or, at most, very weakly first order. This issue is 
generally important, but takes a particularly significant 
role in the context of the local quantum critical point 
(LQCP) The latter has been invoked to explain 

some of the non-Gaussian quantum critical properties ob- 
served in the heavy fermion metals P, |^, 0, Hol . 

In this paper, we study the dynamical competition be- 
tween the RKKY and Kondo interactions using an ex- 
tended dynamical mean field theory (EDMFT) [H [H 
We analyze the global phase diagram of the model 
by carrying out EDMFT studies that cover not only 
the paramagnetic phase but also the antiferromagnet- 
ically ordered phase. Our analysis involves a numeri- 
cal approach to the self-consistent EDMFT equations. 
To understand the nature of the quantum phase tran- 
sition requires a study at sufficiently low temperatures, 
and this is in general very difficult to achieve numeri- 
cally. We have overcome this difficulty by focusing on 
a Kondo lattice model with Ising anisotropy. We are 
able to reach temperatures of the order of about 1% of 
the bare Kondo scale, adapting a quantum Monte Carlo 
(QMC) method|l,4ll(il|. Our most important conclusion 
is that the zero-temperature transition is indeed contin- 



uous and locally critical. We have also determined in 
some detail the dynamics in the ordered and paramag- 
netic phases. 

The model is specified by the following Hamiltonian: 



-^JKS..Se. + ^(/.,/2)5f5|. (1) 



Here, c]g. creates a conduction electron of spin projection 
cr at the i-th site, and Si and s^.i represent the spins of 
the S = ^ local moment and conduction electrons re- 
spectively; tij is the hopping integral [corresponding to 
a band dispersion ] and Jk the Kondo coupling. The 
chemical potential (/i) is chosen so that the system is less 
than half-filled; all the phases we will consider are metal- 
lic. Finally, lij denote the Ising-exchange interactions 
between the local moments. The corresponding Fourier 
transform, /q, is most negative at an antiferromagnetic 
wavevector Q (/q = — /). To access a LQCP, wc will 
consider two-dimensional magnetism |^ and take the 
RKKY -density-of-states in the form: 
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(2) 



where O is the Hcavisidc function. In addition, we con- 
sider a featureless po{e) = ^ ^k)- 

We seek to determine the transition into an antifer- 
romagnetic state. To incorporate an order parameter, 
we separate the inter-site exchange interaction term into 
its Hartree and normal-ordered parts [13, [13: Hi = 
E., [(1/2) : Sf :: 5| : +{SpS, - (l/2)(5f , 
where : Sf := S'f - (S'f). In the EDMFT mapping to 
a self-consistent impurity model, the first term on the 
RHS of Hi leads to a coupling of : S"^ : (at the selected 
site) to a fluctuating magnetic field. Up to constant 
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terms, the effective impurity action takes the form: 

Simp ^ Stop + / dT[hioc S'{t) + JkS{t) ■ Sc(r)] 
Jo 
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where (3 = 1/kBT, Stop describes the Berry phase of the 
local moment, and /iioc, Gq ^, and Xo^ the static and 
dynamical Weiss fields. These fields arc subject to the 
self-consistency condition: 
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with the aid of Dyson-like equations: 'S„{lu) = Gq ;^(cj) — 
g;^, M{uj) = xo'(^) + Here A/((l>) and 

E^(w) are the spin and conduction-electron self-energies, 
respectively; toaf = ('5'^)imp is the staggered magneti- 
zation; xioc(w) and Gioc.(t(^) arc the Fourier transforms 
of {Tr[: S\t) :: 5^(0) :])i„p and -(T,[c,(r)ct(0)])i„p, 
respectively. The lattice spin susceptibility is |l2l |: 



x(q,c^) = l/[/q + Af(a;)] 



(5) 



The self-consistent impurity model, Eq. iPJ, can also 
be written in a Hamiltonian form, 



H 



imp 



hoc S'' + Jk ^-^c+Y^ Ep,a 
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+ 9Y.S^ (4>p + 4>lj) + E -^pVp • (6) 
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It describes a local moment coupled not only to a 
fermionic bath (cpcr), but also to a dissipative bosonic 
bath {<j)p) and a static field. Here Ep„, Wp and g arc such 
that Ep^Wpgyiio'-wD = -Xo\^), Ep^/(^-Epa) - 
Go,(t(w). In order to understand the quantum phase 
transition, we need to reach sufficiently low tempera- 
tures; this is made possible by projecting Eq. © to a 
form that does not involve vastly different energy scales. 
For a Kondo lattice, the ordered staggered magnetization 
should be mostly associated with the localized moments. 
Its main effect in Eq. (jSJ is to generate, in addition to a 
finite hioc, a Zeeman splitting in the fermionic bath. For 
featureless conduction electron density of states we arc 
considering, however, this splitting should yield a negli- 
gible difference between the fermionic density of states at 
the chemical potential for the two spin species. Denoting 



d{LL} — Ep„) = Nq, we can then adopt the same canon- 
ical transformation "Hl^p = WTiim^U within a bosoniza- 

imp 



tion procedure, as used in Ref. 
into 



, to transform H 
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WpCpJcjip , (7) 
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where Ho describes the fermionic bath, is the z- 
component of the fermionic bath spin density, and T and 
Tz are respectively determined by the spin-flip and lon- 
gitudinal components of the Kondo coupling. The cor- 
responding partition function is Z^^^ ~ Tr exp[— ^/j^p] 
where the trace is taken in the spin space and 



SLp = dT[hloc S^t) +rS^-l r dT'S^T)S%T') 

Jo ^ Jo 

y^{Xo\r-T')-K,{T-T'))\. (8) 

Here, ICc{iuJn) = Kc|'-^ti| (i^c ^ ^l^o) a retarded inter- 
action from integrating out the fermionic bath. 

The effective impurity model is solved u sing a quan- 
tum Monte Carlo (QMC) method 0, E [ig^ We it- 
crate on h\oc and Xq ^(t), through the self-consistency 
Eq. Q). In the following wc will consider Kc = t^, and 
A^or ~ 0.19, corresponding to the bare (normalized) 
Kondo scale NqT^ = 0.17. The number of time slices 
in the Trotter decomposition, Nt, was chosen such as to 
make the discretization error small enough. We found 
that Nt > 6 /3 is appropriate. In our simulations TV^ 
varied between 32 for the higher temperatures and 512 
for the lowest one. We used 10^ MC steps per time slice 
and performed between 5 and 25 iterations. 

That antiferromagnetism develops for sufficiently large 
RKKY interactions is seen in Fig. ^a), which shows the 
temperature dependence of the staggered magnetization 
(toaf) as a function of temperature at / = 3T]^ and 
/ 1.5T^. The order parameter toaf becomes non- 
zero at Tn « 0.29rO. and Tn ~ 0.05TO., respectively. 
The numerical data are consistent with a jump in toaf 
at Tn, suggesting that the finite temperature magnetic 
transition is first order. This is further supported by the 
T-dcpendence of the static local susceptibility shown in 
Fig.^b). It is seen that XIoc(T'n) is finite whereas, for an 
RKKY density of states of Eq. (0) , it must be divergent 
if the transition were second order [following arguments 
similar to those given after Eq. ®] 113 ■ 

The dependence of the Ncel temperature on the RKKY 
interaction is shown in Fig.l^^a). It is seen that Tn van- 
ishes at Ic/T'^ « 1.2. The staggered magnetization at the 
lowest temperature (T = O.OllT]^) also vanishes continu- 
ously at the same value Ic as shown in Fig.E^b). Finally, 
the size of the jump of toaf at Tn decreases with the lat- 
ter as the strength of the RKKY interaction is reduced. 
The extrapolated value of the jump also goes to zero at 
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FIG. 1: Temperature dependence of the staggered magnetiza- 
tion mAF (a) and the normalized static local susceptibility (b) 
for the values of the RKKY interaction 1 IT% = 3.0 (squared 
curves) and 1.5 (diamond curves). 




FIG. 2: (a) The Neel temperature (solid dots) and E'{^^ (open 
diamonds) as a function of I /T'^, the ratio of the RKKY in- 
teraction to the bare Kondo scale; (b) The staggered magneti- 
zation (solid dots) as a function of of I/Tg-, at T = O.OllTj^. 
E*^^ characterizes the coherence scale of the paramagnetic 
solution (see the main text). The solid and dashed lines are 
guides to eyes. 

Ic- These observations strongly suggest that the zero- 
temperature transition is second order. 

We also show in Yig.\^E^^^, the coherence scale of the 
Kondo lattice as a function of /. It may be seen that 
E*^^ is finite for I < Ic but decreases with increasing 
/; within our numerical uncertainty, it vanishes precisely 
at Ic- -Ejqj. characterizes the crossover from the low en- 
ergy Fermi liquid regime, in which the local susceptibil- 
ity has a Pauli form, to the high energy quantum critical 
regime where the local susceptibility is logarithmically 
singular. This scale is determined numerically from the 
frequency dependence of the local susceptibility as dis- 
cussed in Ref. 0- In this range (/ < Ic), the solution 
derived from starting the iteration with a zero trial hioc 
is identical to that with a finite trial /iioc- That the val- 
ues of Ic determined from the paramagnetic and ordered 
sides are the same provides an additional evidence for the 
second order nature of the quantum phase transition. 

We wish to stress a subtle issue in the determination 
of Ic coming from the paramagnetic side. If we apply the 
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FIG. 3: Curie constant C in the nominally "paramagnetic" 
solution as a function of I /T%. Inset: the local susceptibility 
vs. Matsubara frequency tOn for I = 2T£- at various temper- 
atures; the jump at c<j„ = measures C/T. 

paramagnetic EDMFT equations (i.e., with hioc — 0) to 
I > Ic, we find that they would converge to a nominally 
"paramagnetic" solution (at least for / < BT]^). These 
solutions are, however, characterized by free residual (but 
not ordered!) moments, as signaled by a C/T divergence 
in Xioc(T') at low temperatures; the constant C is an ef- 
fective Curie constant. This Curie constant manifests 
also in the behavior of xioc('^ri) at a;„ = 0, as illustrated 
(for a particular value / = 2T^ :s> Ic) in the inset to 
Fig- El Xioc(wn = 0) jumps by an amount C/T compared 
to the value extrapolated from xioc(w„) at finite a;„ [l8| . 
(Without going to low temperatures, this jump would 
have been hard to see.) From the dependence of C on / 
(Fig. 131 main plot), the Curie constant is seen to vanish 
as / approaches Ic from above. The finite Curie constant 
for all I > Ic shows that, in this region, the nominally 
"paramagnetic" solution is unstable and the true ground 
state is the antiferromagnetic solution. 

That E*^^ has to vanish at Ic can in fact be under- 
stood once the continuous onset of Tn and ttiaf at Ic is 
established. A second order magnetic transition at zero 
temperature means that xi^,^ = 0,T = 0) diverges at 
Ic- Using 

X(Q, LO^0,T^0) = 1/[M{lu = 0, T = 0) - /] (9) 

[which can be seen by inserting Iq = —I into Eq. (Q], 
this imphes that M{iu = 0, T = 0) = /. It then fol- 
lows from the self-consistency Eq. (j4b|l . together with 
Eq. ISJ, that the local susceptibility xioc('^ = 0, T = 0) 
also diverges 0| at I = Ic- This places the local Kondo 
problem, as specified by Eq. © with h\oc ~ 0, to be at 
its own critical point |l7| . In other words, the criticality 
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FIG. 4: The normalized local susceptibility vs. Matsubara 
frequency, at T = O.OUT]^. The legend specifies //T£,. 

of the Kondo physics is embedded in the criticahty of 
the magnetic transition in the lattice; the corresponding 
coherence energy scale, E*^^, vanishes. 

The overall phase diagram is now specified by Fig.l^a). 
While the first order nature of the finite temperature 
transition is an artifact of EDMFT, reflecting the fact 
that EDMFT contains no spatial anomalous dimension, 
the LQCP at zero temperature is expected to be ro- 
bust 0, ^3- establishing the continuous nature of 
the zero-temperature transition, we can identify a finitc- 
temperaturc quantum-critical "fan" that is anchored by 
the (T = 0) LQCP at / = /c, even in the absence of 
a correct understanding of the finite temperature transi- 
tion. This quantum critical fan can be used to describe 
the non-Gaussian quantum critical behavior observed in 
heavy fermion metals. 

To better understand the dynamics in the ordered 
phase, we have also calculated the frequency dependence 
of the local susceptibility. Fig. ^ plots xiod^^n) vs cun 
for various values of //T^, at T = O.OllT]^. On the 
ordered side, the peak susceptibility, x(Q, w„), saturates 
to a finite value as a result of the finite magnetic or- 
der parameter. Correspondingly, the local susceptibility, 
Xioc('^n), which is equal to the q-average of x(q, also 
saturates. The saturated value of xioc(^n) increases as 
the RKKY interaction and, correspondingly, /iioc, is re- 
duced. At I = Ic ~ 1.2T]f, hioc goes to zero and the 
local susceptibility becomes logarithmically singular, the 
form that is associated with the critical point of this self- 
consistent impurity model; the result at I = Ic agrees 
with that of the analysis from the paramagnetic side. 
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